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ABSTRACT 

From two very simple axioms: (1) that AGN activity traces spheroid formation, and 
(2) that the cosmic star-formation history is dominated by spheroid formation at high 
redshift, we derive simple expressions for the star-formation histories of spheroids and 
discs, and their implied metal enrichment histories. 

Adopting a Baldry-Glazcbrook initial mass function wc use these relations and 
apply PEGASE.2 to predict the z = cosmic spectral energy distributions (CSEDs) of 
spheroids and discs. The model predictions compare favourably to the dust-corrected 
CSED recently reported by the Galaxy And Mass Assembly (GAMA) team from the 
FUV through to the K band. The model also provides a reasonable fit to the total 
stellar mass contained within spheroid and disc structures as recently reported by 
the Millennium Galaxy Catalogue team. Three interesting inferences can be made 
following our axioms: (1) there is a transition redshift at z « 1.7 at which point 
the Universe switches from what we refer to as "hot mode evolution" (i.e., spheroid 
formation/growth via mergers and/or collapse) to what we term "cold mode evolution" 
(i.e., disc formation/growth via gas infall and minor mergers); (2) there is little or no 
need for any pre-enrichment prior to the main phase of star-formation; (3) in the 
present Universe mass-loss is fairly evenly balanced with star-formation holding the 
integrated stellar mass density close to a constant value. 

The model provides a simple prediction of the energy output from 
spheroid and disc projenitors, the build-up of spheroid and disc mass, and 
the mean metallicity enrichment of the Universe. 

Key words: galaxies: formation — galaxies: evolution — galaxies: bulges — galaxies: 
ellipticals — galaxies: star formation — galaxies: spiral 



1 INTRODUCTION 

The vast majority of galaxies can be adequately described 
as consisting of a compact smooth spheroidal component 
containing a predominantly pressure-supported old [a/Fe]- 
enhanced stellar population, and/or an extended flattened 
star-forming disc component containing intermediate and 
young stars with a wide range in metallicities, having smooth 
rotation and embedded in an extensive gaseous cold gas 
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disc. Exceptions exist, most notably the dwarf populations 
which, while dominant in terms of number density, actu- 
ally contribute only a modest amount to the baryon budget 
at the present time (< 16 per cent; Driver 1999; Geller et 
al. 2012). This dichotomy of galaxies into spheroids and discs 
has been known for over a hundred years stretching back to 
even before the confirmation that galaxies are external sys- 
tems (e.g., Hubble 1926, 1936; Zwicky 1957 and references 
therein). To some extent this dichotomy has been recently 
"rediscovered" , through the statistical studies of large pop- 
ulations as a galaxy bimodality (Strateva et al. 2001; Baldry 
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et al. 2004). Driver et al. (2006) argued that this bimodality 
is better interpreted in terms of the earlier bulge-disc di- 
chotomy and advocated routine structural decomposition as 
vital (e.g., Allen et al. 2006; Simard et al. 2011; Lackner & 
Gunn 2012) to directly trace the independent evolutionary 
histories of the spheroidal and disc components. 

Numerical models of galaxy formation struggle to pro- 
duce realistic galaxy systems with a tendency to form overly 
cuspy cores and difficulty in maintaining extended disc 
structures with a high axis ratio (White & Navarro 1993; 
Navarro & Steinmetz 2000; Abadi et al. 2003; House et 
al. 2011). Both problems are likely to be connected to the 
different fundamental properties of the dark matter and the 
baryons, and in particular their ability to experience pres- 
sure and their ability to dissipate energy. In the core regions 
the gravitational coupling of the baryons with the dark- 
matter may allow it to exhibit a pseudo-pressure, whereas 
in the outer-regions the ability of the baryons to dissi- 
pate energy on a timescale which is faster than the free- 
fall timescale may allow for the formation of a thin rotating 
baryonic disk. This picture while simple to articulate has 
proven extremely hard to simulate, with the need to parti- 
tion and redistribute the angular momentum in a quite spe- 
cific manner to result in galaxies with realistic appearances. 
In particular merger events are extremely disruptive to this 
process, imparting both energy and angular momentum to 
the baryonic disc, which is easily disrupted or 'plumped 
up' (see Barnes & Hernquist 1992 for extensive discussion 
and early references on this topic, also Hopkins et al. 2009 
for updated simulations on the survivability of discs during 
merger events). In general the greater the merger-rate the 
more bulge-dominated the final gala^cy population appears. 

More contemporary hydrodynamical simulations (e.g., 
Governato et al. 2010; Agertz, Romain & Moore 2011; Scan- 
napieco et al. 2011; Domenech-Moral et al. 2012) are now 
starting to show significant success at producing realistic 
"looking" bulge-disc systems by incorporating a greater level 
of cold gas infall than previously assumed, as argued ear- 
lier by Keres et al. (2005) and Dekel et al. (2009). These 
focused hydrodynamical studies, however, are inevitably ex- 
tracted from numerical simulations with particularly qui- 
escent merger histories, suggesting such systems should be 
the exception rather than the norm. Hence, while numerical 
simulations of the development of the dark matter haloes 
find a continual process of halo merging, it appears that the 
baryons and what we identify as galaxies (baryonic conden- 
sates), might not develop in the same way. Martig & Bour- 
naud (2009) argued that feedback from low and intermedi- 
ate mass stars can contribute significantly to the redistribu- 
tion of mass from the bulge to the disc through extensive 
(or even excessive) mass-loss. This baryonic outflow could 
help alleviate the problem of excessive bulge-formation by 
allowing some fraction of the collapsed stellar mass (up to 
~ 50 per cent), to return to the halo and contribute to the 
later growth of a disc — thereby coupling bulge and disc 
growth. However this mechanism also relies on a fairly qui- 
escent merger history during later times, and does not easily 
explain the broad diversity of bulge-disc ratios seen. As an 
aside, simulations have also demonstrated that baryonic out- 
flows from the core regions can help provide a plausible ex- 
planation to the core-cusp problem (Governato et al. 2012; 
Zolotov et al. 2012). Clearly feedback and infall are both 



crucial processes, whose motivation is as much driven by 
the requirement to produce realistic looking images as by 
fundamental physics, and which are both more effective if 
the merger rate is either low or at least confined to earlier 
epochs. 

As argued in the opening paragraph, we advocate a 
more heuristic approach where we put aside the issue of 
dark-matter assembly and start by asking whether the di- 
chotomy of galaxy structure is best explained by two distinct 
formation mechanisms. Following the earlier discussion and 
lessons learnt from the simulations, the obvious two mech- 
anisms can loosely be termed as a hot and cold mode. In 
the hot mode spheroids are formed early and rapidly via 
dynamically hot (turbulent) processes (collapse, fragmenta- 
tion, and merging). In the cold mode discs are formed more 
slowly, from an extended quiescent phase of cold gas infall 
regulated by internal feedback (i.e., supernova). This basic 
concept is of course not new (e.g., Larson 1976; Tinsley & 
Larson 1978) but has laid dormant for sometime overshad- 
owed by the dominance of merger-driven evolution. However 
the revival is also being championed via a series of semi- 
analytic studies by Cook et al. (2009; 2010a; 2010b, see also 
Dekel et al. 2009) inspired by behaviour seen in numerical 
simulations in which an initial rapid hot merger phase is typ- 
ically followed by a more quiescent phase of accretion (see 
also L'Huillier, Combes & Semelin 2012). 

The two-phase model is both obvious (given the bulge- 
disc nature of galaxies) and controversial, as it marginalises 
the merger rate required for dark matter assembly to earlier 
epochs than simulations typically suggest. This low merger 
rate is arguably corroborated by the local studies of dynami- 
cally close pairs (in particular see Patton et al. 2002; De Pro- 
pris et al. 2005, 2007, 2010) — although a correct derivation 
of the merger rates requires a robust understanding of the 
merger timescales, which are currently poorly constrained. 
Perhaps more compelling, however, is the result that only 
40% of the present day stellar mass resides in spheroidal sys- 
tems (Driver et al. 2008; Gadotti 2009; Tasca & White 2011). 
A key inference is then: // discs are destroyed/thickened 
during mergers, yet the majority of stellar mass resides in 
discs, the dominant formation mechanism cannot be merger- 
driven, but presumably the more quiescent process of cold 
gas accretion. This statement becomes more profound when 
one realises that the stellar mass in discs today only mea- 
sures that unaffected by mergers, and that some of the stars 
currently in spheroidal systems may have originally formed 
within discs via cold accretion prior to a merger event. In 
some bulge formation scenarios, star-formation via merging 
is dispensed with altogether and replaced by the migration of 
massive star-formation clumps formed within deeply turbu- 
lent discs (e.g., Elmegreen, Bournaud & Elmegreen 2008). 
This potentially relegates the stellar-mass build-up driven 
by mergers to be in the — 40 per cent range by mass. 
Clearly mergers do occur at all redshifts and similarly discs 
may form, be disrupted, and reform at any redshift. Re- 
cently L'Huillier, Combes & Semelin (2012) reported that 
77 per cent of a galaxy's mass is formed via gas accretion 
and 23 per cent via direct merging from simulations. Other 
empirical studies also seem to suggest that the bulk (~ 70 
per cent) of the stellar mass is mostly assembled by z ~ 1, 
again marginalising the role of late time major mergers (e.g., 
Bundy et al. 2004; Brown et al. 2007, 2008). 



Two-phase galaxy evolution 3 



Focused studies of nearby galaxies are also unveiling 
significant levels of gas accretion in some nearby systems 
(Sancisi et al. 2008) and studies of the very rapid evolu- 
tion of galaxy sizes have argued (e.g., Graham et al. 2011) 
that the compact elliptical systems seen at intermediate red- 
shift (1.4 < z < 2.5) by Daddi et al. (2005) and Trujillo et 
al. (2006) (see also Bruce et al. 2012) might represent the 
naked bulges of present day spiral systems. 

In essence the two-phase model is an attempt to high- 
light, conceptually, the possibility of a distinct change in 
the primary galaxy formation mechanism occurring at some 
transition redshift from an era where the dominant mode 
is major mergers leading to spheroid formation, to an era 
where the dominant mode is accretion leading to disc for- 
mation. 

At early cosmic epochs we see a prevalence of distinct 
phenomena, in particular highly asymmetrical morphology 
in massive/luminous systems (Driver et al. 1998; Conselice, 
Blackburn & Papovich. 2005; Ravindranath et al. 2006) and 
significantly increased AGN activity (Fan et al. 2001,2003; 
Groom et al. 2004; Richards et al. 2006). AGN activity is di- 
rectly linked to the formation and growth of the associated 
super-massive black holes (SMBH; Hopkins et al. 2008a) 
which in turn is linked to spheroid formation via the well 
established SMBH-bulge relations (see, for example, the re- 
view by Ferrarese & Ford 2005 or the recent near-IR SMBH- 
bulge luminosity relation by Vika et al. 2012). Recent stud- 
ies also argue, from more direct empirical evidence, that 
AGN activity is almost always coincident with massive star- 
formation and that the two-processes do indeed appear to 
occur hand-in-hand (e.g., Rafferty et al. 2011). This AGN- 
SMBH-bulge connection therefore implies a clear timescale 
for the formation of the spheroid systems (see also Pereira 
& Miranda 2011 for a similar argument, albeit applied in 
the opposite direction). 

In Section 2 we describe the z — empirical data de- 
scribing the cosmic spectral energy distribution (GSED) of 
spheroids and discs as recently reported by the Galaxy And 
Mass Assembly team (Driver et al. 2011, 2012). In Section 3 
we take the above arguments to their natural conclusion and 
use the AGN-SMBH-Bulge connection to define the indepen- 
dent star-formation history of spheroids and assign the resid- 
ual star-formation, implied by the cosmic star-formation his- 
tory, to describe that of the discs. In Section 4 we use our 
star-formation histories to produce predictions of the CSED 
of spheroids and discs, and in Section 5 compare the predic- 
tions to the data. 

Throughout we use Hq—YO h km s^^ Mpc^'^ and adopt 
Qm = 0.27 and = 0.73 (Komatsu et al. 2011). 



2 THE Z ^0 COSMIC SPECTRAL ENERGY 
DISTRIBUTION 

In Driver et al. (2012) we reported the empirical measure- 
ment of the cosmic spectral energy distribution (GSED) in 
the nearby Universe, corrected for dust attenuation, and 
spanning the wavelength range from 0.1 to 2.1 /^m, i.e., the 
regime over which direct stellar-light dominates. These data 
were derived from the combination of the GAMA spectro- 
scopic survey, currently underway on the Anglo- Australian 
Telescope (Driver et al. 2011), coupled with reprocessed and 



aperture matched data from GALEX, SDSS, and UKIRT 
LAS (see Seibert et al., in prep, and Hill et al. 2010). Driver 
et al. (2012) also provided the GSED subdivided according 
to spheroid-dominated and disc-dominated systems. The di- 
vision into spheroid and disc dominated was achieved via 
visual classification, as neither a simple colour nor Sersic in- 
dex division appears to cleanly separate the two populations 
(see also Kelvin et al., 2012, figure 20). 

The sample originated from a common volume of 2.8 x 
10^ (Mpc/h)^ over the redshift range 0.013 < z < 0.1. Al- 
though the GAMA survey currently contains about 180,000 
galaxies with known redshifts, the adopted redshift range 
significantly reduces the sample size to around 10,000. It 
also simplifies and removes any luminosity bias arising from 
large scale clustering as the sample is pseudo- volume limited 
around the L* region — i.e., those galaxies which dominate 
the luminosity density measurements. 

As the GAMA survey lies entirely within the Sloan Dig- 
ital Sky Survey, the GAMA CSED was renormalised to the 
full SDSS survey area. This reduces the cosmic/sample vari- 
ance from around 15 per cent to 5 per cent (using the formula 
for estimating cosmic variance given by eqn. 4 of Driver & 
Robotham 2010). 

2.1 The CSED of spheroids and discs 

The final CSED values we use here are the spheroid- 
dominated and attenuation corrected disc-dominated values 
taken directly from Table 7 of Driver et al. (2012). 

As dust attenuation is such a crucial issue it is worth 
mentioning the genesis of the corrections used by the GAMA 
team. First, the dust correction is only applied to the disc- 
dominated data and the spheroid population is assumed 
dust free (e.g., Rowlands et al. 2012). Second, the correc- 
tions are based on the radiative transfer models of Tuff's 
et al. (2002) and Popescu et al. (2011) which have been 
fine-tuned to the multiwavelength (FUV-far-IR) data of 
NGG891, and incorporate three distinct dust components; 
an extended low opacity disc, a compact high opacity disc 
and dust clumps. This fiducial model is calibrated to the 
galaxy population at large by modifying the B-band central 
face on opacity until the predicted variation of flux with in- 
clination matches the trend of M* with inclination seen in 
the Millennium Galaxy Catalogue data (Driver et al. 2007). 
This calibrated model was then used to derive the combined 
face-on and inclination dependent correction for a popula- 
tion of galaxies averaged over all viewing angles and over 
a wavelength range of 0.1 — 2.1 /im (Driver et al. 2008). 
This photon escape fraction (varying from 24 per cent in 
the FUV to 89 per cent in the Jf-band) was then used to 
correct the GSED of disc-dominated systems to the intrinsic 
GSED, which we use here. 

The CSED of spheroid-dominated and disc-dominated 
galaxies, however, is not precisely what we require, as some 
proportion of the CSED flux in the disc-dominated class 
may be coming from the central bulges. Likewise, some pro- 
portion of the flux in the spheroid-dominated class may be 
due to faint discs. In order to assess how much of a problem 
this might be, we can compare the ratio of the _ftr-band lumi- 
nosity densities of the spheroid-dominated to non-spheroid 
dominated samples, to the ratio of the stellar mass densi- 
ties of bulge-|-elliptical systems to disc systems from Driver 
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et al. (20070. This test assumes that the K-haiid lumi- 
nosity is a suitable single-band proxy for stellar mass. We 
find reasonable agreement (within 12 per cent), suggesting 
that a comparable amount of flux needs to be redistributed 
in either direction. In detail, the if-band luminosity densi- 
ties are 1.2 and 2.2 (xlO^" h W Mpc"^) for the spheroid- 
dominated and non-spheroid-dominated populations respec- 
tively (taken from Table 7 of Driver et al. 2012). Meanwhile 
the stellar mass densities for spheroids (bulge-(-ellipticals) 
and discs are 2.9 and 4.7 (xlO*h Mq Mpc~^) respectively 
(taken from Table 1 of Driver et al. 2007) . This gives agree- 
ment to ~ 10 per cent and suggests that for the moment we 
can adopt the following approximation: 

elliptical-l-bulge CSED f» Spheroid-dominated CSED 

disc CSED f» non-Spheroid dominated CSED 

In due course all galaxies at z < 0.1 in the GAM A 
survey will be decomposed into bulge and disc components 
to enable a direct derivation of the true spheroid and disc 
CSEDs. 



3 THE STAR-FORMATION HISTORY OF 
SPHEROIDS AND DISCS 

The local CSED should be a predictable quantity if the cos- 
mic star-formation history (CSFH) is known, the initial mass 
function (IMF) is universal and known, and a plausible stel- 
lar evolution code applied. Of course this is not quite so sim- 
ple and in particular the metallicity enrichment adopted will 
significantly modify the predicted CSED shape. Upcoming 
papers will explore these issues in more detail, but here we 
wish to construct a basic first-look model and focus on the 
viability of the hypothesis that galaxy formation progressed 
in two fairly distinct phases: rapid spheroid formation fol- 
lowed by more quiescent disc growth. 

In order to construct a model of the present day 
spheroid and disc CSEDs we need not just the CSFH, but 
the CSFH sub-divided into spheroids and discs. These CSED 
predictions can then be compared to the data from Section 
2. 

The existence of various super-massive black-hole bulge 
relations (e.g., Ferrarese & Ford 2005), provides the obvious 
smoking gun, as it couples SMBH growth to bulge growth. 
This is because SMBHs are believed to grow via mergers, 
resulting in an active-galactic nucleus phase (Hopkins et 
al. 2006). The growth of spheroids is therefore, arguably, 
mirrored via the more readily observable AGN activity his- 
tory. This logical connection, from a correlation to causality, 
is the key assumption underpinning our model and forms the 
first of our two axioms. In the recent study by Richards et 
al. (2006), the integrated AGN activity versus redshift was 
reported and, ignoring any significant lag (in either spheroid 
formation or AGN activity), can be used as a proxy in shape 
for the spheroid cosmic star-formation history. 

The amplitude of the spheroid SFH can be set from 

^ Tlie Driver et al. (2007) study is based on bulge-disc decomposi- 
tions of the Millennium Galaxy Catalogue data (Liske et al. 2003; 
Driver et al. 2005) described in full in Allen et al. (2006) 



comparison of the AGN activity shape to the global CSFH. 
For our second axiom we elect to maximise the spheroid 
CSFH by setting the amplitude as high as possible without 
exceeding the global CSFH (i.e., a maximal spheroid forma- 
tion scenario). Conceptually then, the heart of the two- phase 
model can be defined empirically from two axioms: 

1) AGN activity traces spheroid growth 

2) Spheroid formation dominates at high-z 

As the above two axioms are already constrained by em- 
pirical data, this provides a zero-parameter starting point for 
the two-phase model — bypassing the need for any initial 
conditions or detailed numerical simulations. Fig. [1] (upper) 
shows the fit (solid curve) to the cosmic star-formation his- 
tory data (grey data points) taken from Hopkins & Bea- 
com (2006; see their figure 2a and table 1 column 2). This 
adopts the parametric form defined by Cole et al. (2001) 
and where the UV data has been calibrated to a modified- 
Salpeter (1955) IMF. The data describing the AGN luminos- 
ity density are taken from Richards et al. (2006) and rescaled 
such that the peak of the AGN luminosity density lies on 
the CSFH curve (requiring an arbitrary multiplication by a 
factor of 3.51 x lO^M0/yr~^LQ). Immediately noticeable is 
the apparent discrepancy/uncertainty at very high redshift. 
Particularly as the axioms above require that at the very 
highest redshift the CSFH and AGN activity curves should 
have the same form. To some extent the evident discrepancy 
simply reflects data uncertainty, as the dust corrections on 
the CSFH at high-z are poorly constrained with significant 
ongoing debate as to the true shape of the CSFH at the 
highest redshifts. For example, measurements of the star- 
formation history based on gamma-ray bursts (Yuksel et 
al. 2008; Kistler et al. 2009) often find higher star-formation 
rates. Likewise the incidence of dust-obscured AGN is an 
equally hotly debated topic (e.g., PoUetta et al. 2006; Tri- 
ester et al. 2010). Most recently Behroozi, Wechsler 
& Conroy (2012) argue that the compendium by 
Hopkins & Beacom (2006) potentially leads to an 
over-estimate of the cosmic star-formation history 
at very high redshifts as the pre-2006 UV luminos- 
ity densities may have been over-estimated and find 
a modified CSFH which agrees very well with the 
high-z AGN data (see their figure 2). 

When presented as a function of time (Fig.[T] lower), it 
is clear that this discrepancy is not actually that significant, 
and the accumulated stellar mass able to be formed during 
this high-z interval is small compared to subsequent mass 
growth. We can now trivially fit the modified-Richards data 
to derive the star-formation history of spheroids. This fit 
can then be subtract from the global star-formation history 
and in turn fit to recover the implied disc formation history. 
The resulting expressions are given below and represent the 
star-formation rate (p) versus time in Gyrs (tcyrs) since the 
Big Bang for the spheroidal (S) and disc (D) populations: 



PS = Cl-03 X 10"^/ii].7( 



tGyraho.r 



v8.57 / 21.86 , , , 

) exp(-- — ) (1) 

EGyrsfto.T 



pz, =ei.80 X lO-^hU-^^r^\M-T^^) (2) 

tGyrs/lO.7 IGyrsAlO.T 

where ^ is the IMF multiplier as given in Table[T] The guide- 
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Figure 1. [upper panel) the cosmic star-formation history from Hopkins & Beacom (2006), Table 2 (black line) and the actual data (grey 
points) calibrated to a modified-Salpeter IMF (see Table [T]|. Overlain are the QSO luminosity density data from Richards et al. (2006), 
Figure 20. The QSO luminosity density is scaled until the peak aligns with the peak of the CSFH. {lower panel) the same data but now 
shown with the ordinate in units of age. Overlain are our parametric fits to these data which provide our inferred CSFHs for the Spheroid 
and Disc systems. On both panels we also show the SFR derived from the SDSS/GALEX FUV luminosity function given in Robotham 
& Driver (2006) converted to a modified-Salpeter A IMF (see Table [T] or Hopkins & Beacom 2006). 



line error at any particular time should be taken as --^ ±25% 
based on the scatter of the original data shown on Fig. [1] 
(lower) i.e., ~ 70%of the grey data points lie on the grey 
shaded regions. At this point it is also worth highlighting 
that while the AGN data (mauve points) place a strong con- 
strain on the shape of the spheroid star-formation history 
the normalisation is very uncertain and based on the rela- 
tively limited high-z cosmic star-formation history data and 



as we shall see in Section. 4 will require a modest downward 
correction of 25% to match the observed CSED (i.e., within 
the grey shaded error bounds). 

As an alternative for the spheroid population, one could 
instead use the total CSFH for t < 3Gyrs combined with the 
AGN luminosity density data for t > 3Gyrs, which is given 
by: 
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Table 1. CSFH multiplication factors {(,) for various IMFs. 
IMF Multiplier 

i 

Salpoter (1955) xO 
modified-Salpeter (Hopkins & Beacom 2006) xl.O 
Baldry & Glazebrook (2003) xO.7 
Kroupa (1993) xl.7 
Kroupa (2001) xO.85 
Chabricr (2003) xO.85 

1 « QO 1 « CO 



PS2 = C2.72 X lO^^hU-^^^r' eM-j^^) (3) 

lGyrsfto.7 tGyrs'lO.7 

As stated these CSFHs are calibrated, via the UV data, 
to the modified-Salpeter IMF used by Hopkins & Beacom 
(2006) which was first laid down as Sal A by Baldry & Glaze- 
brook (2003). To convert to other IMFs one needs to multiply 
by the factor (^) shown in Table [T] The simple expressions 
above are shown in Fig. [1] (lower) in red (spheroid), blue 
(disc), and green (spheroid -I- disc) and provide a good fit to 
the data given the accuracy to which the data are known. 
These equations now provide a blueprint for the formation 
of the present day spheroids and discs over the full age of the 
Universe, leading to a clear prediction of the stellar energy 
output and stellar mass growth. 

Of particular interest should be the transition point 
around 4.2Gyrs [z ~ 1.6) from which point star-formation 
resulting in disc growth dominates over star-formation re- 
sulting in spheroid growth. This suggests a key epoch at 
which the Universe switches from merger dominated evolu- 
tion to accretion dominated evolution and meshes very well 
with the evident change in the morphological appearance 
of galaxies from highly disturbed to more ordered systems 
at z ~ 1.5 (see Driver et al. 1998; and also van den Bergh 
2002). 



4 CONSTRUCTING THE MODEL 

With the CSFH of spheroids and discs defined, we now build 
our empirical two-phase model adopting "vanilla" choices at 
every opportunity. In summary the key inputs and assump- 
tions to the model are: 

1) The star-formation history shown in Fig. 1 as discussed 
in section 3. 

2) The adoption of a Universal IMF, in this case Baldry & 
Glazebrook (2003, henceforth BG03). 

3) The adoption of pegase.2 (Fioc & Rocca- 
Volmerange 1997; 1999) to model the spectral output 
of the evolving stellar population (using default options 
throughout). 

4) The assumption that the gas-phase metallicity increases 
linearly with star-formation from Z — 0.0 to Z = 0.030 for 
spheroids and to Z = 0.010 for discs, with no time lag (i.e., 
instantaneous enrichment). 



4.1 Metal/chemical enrichment history 

Perhaps the most uncertain of the above list is the ap- 
propriate metallicity history to adopt. Here we have been 
guided by the study of Tremonti et al. (2004) to 
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Figure 2. The adopted (gas-phase) metallicity for spheroids 
(red solid line) and discs (cyan solid line) as a function of time. 
Solid data points show the accepted mean 2 = values taken from 
Tremonti et al. (2004). Also shown are approximate data values 
from Erb et al., (2006) and Zahid et al., (2011). Note that all 
data have an arbitrary error of ±10 %. The gas-phase metallic- 
ity assumes no lag between star-formation and enrichment. The 
dotted lines show the implied metallicity of the stellar popula- 
tions for the spheroids and discs. The grey line shows the mean 
by mass of the integrated stellar metallicity and the instantaneous 
integrated gas-phase metallicity. 

define our gas-phase metallicity at redshift zero 
to be Z = 0.030 and Z = 0.010 for the spheroids 
and discs. These values were determined by noting 
the metallicity at 10^^ M0 (predominantly spheroids) 
and at IO^Mq (predominantly discs). To convert the 
given 12 + logjg(0///) values to those shown on Figj2] 
we adopt a solar metallicity of Zq = 0.019 with 
12 + \og^Q(0/ H)q = 8.9. We then argue that in the ab- 
sence of other factors the mean metallicity will rise approx- 
imately linearly with the cumulative cosmic star-formation 
history normalised to the present day values. This ignores 
the prospect of either pre-enrichment via, for example. Pop- 
ulation HI stars, or any lag between the star-formation and 
the increase in metallicity (i.e., instantaneous enrichment). 

Conceptually these are loosely consistent with a closed- 
box model for spheroids (i.e., rising to a metallicity close 
to typical yields), and an infall model for discs or one in 
which the disc is gradually growing from a large gas "reser- 
voir" (e.g., perhaps analogous to the "equilibrium model" 
put forward by Dave, Finlator & Oppenheimer 2012). To ex- 
plore the bounds and importance of this metal enrichment, 
however, we also show the CSED predictions using our sim- 
ple evolving metallicity history and for constant metallicity 
at the highest and lowest values. Fig. [2] shows the implied 
metallicity histories derived from Eqns. 1&2, for the two 
populations (as indicated by the red and cyan solid lines). 
Note that the grey lines on Fig. [2] show the com- 
bined metallicity of all stars formed (grey dotted 
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Figure 3. The zero-parameter output, assuming a BG03 IMF for various mctallicity histories as indicated and adopting the star- 
formation histories given by Eqns. 1 & 2. Note that the star-formation rates have be multiplied by a factor of 0.55 to convert from a 
Salpeter (1955) IMF to that of BG03. The data points are transcribed directly from the CSED reported in Driver et al.( 2012) where the 
red points represent spheroid-dominated, the blue disc-dominated and the black the sum of the two. The model lines are for spheroids 
(orange), discs (cyan) and the sum (black). 



line) and the integrated gas-phase metallicity (grey 
solid line). This shows interesting behaviour in that 
the average gas-phase metallicity (of the gas about 
to form stars) peaked at z ~ 2. 

Constraints from the literature on the mean metallic- 
ity at intermediate and high redshift are minimal, however 
we note that Zahid, Kewley & Bresolin (2011) find from a 
sample of 1350 galaxies drawn from DEEP2 that massive 



systems have comparable gas-phase metallicity ai z — 0.8 
to local systems, while low mass systems have a gas-phase 
metallicity reduced by O.lSdex. However Erb et al., (2006) 
find that the implied gas-phase metallicity, for massive, 
i.e., spheroidal-like systems, at z ~ 2 is approximately 
half that at 2; = 0. Both of these results are crudely consis- 
tent with our inferred metallicity history if one equates (as 
we explicitly do), the massive systems to spheroids and low 
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mass systems to discs. Note that one natural byprod- 
uct of this is that as intermediate mass systems have 
both a spheroid and a disc component their systemic 
metallicity will lie somewhere between the two ex- 
tremes and exhibit strong radial gradients as one 
moves from the central spheroid component to the 
outer disc component. As the mean bulge-to-total 
ratio increases fairly smoothly with stellar mass this 
naturally gives rise to the mass-metallicity relation 
(Tremonti et al.2004). Note that the ±10% error ranges 
shown on Fig. [2] are purely indicative as the actual ranges 
are poorly constrained. 



4.2 Stellar population synthesis 

To construct the redshift zero CSED the pegase.2 code 
(Fioc & Rocca-Volmerange 1997; 1999) was used to produce 
a series of single-stellar population (SSP) templates with 
an appropriate range of metallicities (Z = 0.000 to 0.025 in 
0.001 intervals) and with the pegase.2 default steps in ages, 
(i.e., roughly logarithmic from 0-20Gyrs). For all SSP tem- 
plates the star formation is set to a short continuous burst 
over a IMyr period with constant metallicity (leading to the 
formation of 2.0 x 10~^ M© in pegase. 2-norma/iserf stellar 
mass units). These SSP spectra were then combined to cre- 
ate a library of 1 Gyr time averaged spectra from 0-1 Gyr to 
13-14 Gyr in 1 Gyr intervals, and for each metallicity class. 
Note that the 0-1 Gyr bin which dominates the FUV and 
NUV region is extremely hard to model correctly because of 
the rapidly changing UV flux and requires more care. Here 
we take the rather simplistic approach of combining all the 
spectra provided by pegase.2 in the 0-1 Gyr range in the 
following manner i.e., 



1 + 2 + 3+.. . + 10 
10 



+ 20 + 30+. .. + 100 



-+200+300+. 



■ Myr 



- IGyr : 



To create the CSED at any redshift we then sum all pre- 
viously formed populations, aged appropriately, drawn from 
the appropriate metallicity class, and scaled by the required 
star-formation rate. The modelling approach we adopt is 
therefore relatively simplistic and effectively assumes all val- 
ues (star-formation rate, metallicity etc) are held constant 
over a 1 Gyr time period. At this stage we feel this is suf- 
ficient time resolution given the inherent uncertainties in 
the initial assumptions (i.e., the input CSFH and AGN ac- 
tivity data). CSEDs were then derived at all 13 time steps 
and combined to produce simple evolution movies available 
from: 



|http://star- www.st-and.ac.uk/~spd3/modell.gif 
ing metallicity 



evolv- 



[http: / / star-www.st-and.ac.uk/~spd3 / model2.gif 

stant high metallicity 



http://star-www.st-and.ac.uk/~spd3/model3.gif 




stant zero metallicity 



Figure 4. The zero-parameter output for alternative IMFs using 
the evolving metallieity shown in Fig[2] and adopting the star- 
formation histories given by Eqns. 1 & 2. Note that in generating 
the models we modify the input star-formation by factors of X 1.0, 
Xl.3 and xO.7 respectively. 



4.3 Normalisation of the CSEDs 

In order to determine the correct normalisation we need 
to multiply the output pegase.2 SSP spectra which are in 
units of ergs~^ by j^pQ^f^ - Here the factor 10^ scales to 
IGyr bins, the factor 10 converts ergs~^ to W, the wave- 
length is in Angstroms, and the factor 0.002 scales the spec- 
tra to 1 solar mass. In applying Eqns. 1-3 we set ^ to 0.7 to 
correct the CSFHs to the BG03 IMF. 

Finally we allow the normalisations to float by ±25% 
to account for the uncertainty highlighted in Fig. [T] by the 
grey shading, along with uncertainties in the multiplication 
factors in Table [Jl the cosmic variance in the GAM A CSED 
data, and the impact of metallicity on star-formation rates. 
These 1 values are shown in brackets in Figs. [3] & U) 



5 MODELS V DATA 

5.1 The CSED and adopted metallicity 

Fig. [3] shows the direct comparison of our z = Q CSED mod- 
els against the recent GAMA data (Driver et al. 2012), for 
the three assumed metallicity histories (as indicated). The 
top panel, which adopts the evolving metallicity, shows a re- 
markable agreement across the full wavelength range and for 
both the spheroid and disc systems. Note that in achieving 
these fits the spheroid data have been renormalised down- 
wards by 25% which is within the specified range of uncer- 
tainty. The central and lower panels of Fig. [3] show the same 
models except for a constant high or low metallicity. This 
has a negligible impact on the disc-CSED, suggest- 
ing very little dependency on the assumed metal- 
licity evolution for discs (perhaps in part due to 
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Figure 5. The implied build-up of stellar mass in spheroids and 
discs versus recent measurements from the Millennium Galaxy 
Catalogue. Shown in grey are the compendium of data from 
Wilkins, Trentham & Hopkins. 2008a. 



the low value adopted). Conversely the impact on the 
spheroid CSED is quite marked with the CSED tilting ei- 
ther redward or blueward for constant high or constant low 
metallicity respectively. This perhaps lends some argument 
against any very strong pre-enrichment phase as interme- 
diate and low metallicity stars in spheroids are required to 
produce a plausible CSED. The obvious caveat is whether 
the shape of the currently adopted spheroid CSED is signif- 
icantly modified/contaminated by young disc light. 



5.3 Stellar mass history 

Fig.[5]shows the implied build-up of stellar mass in spheroids 
and discs (and combined), as indicated. Note that we show 
both the total cumulative stellar mass formed (dashed lines) , 
along with that remaining based on default pegase.2 as- 
sumptions as to mass-loss (solid lines). Also shown are the 
direct empirical stellar mass measurements from the Mil- 
lennium Galaxy Catalogue (Liske et al. 2003; Driver et 
al. 2005), which includes corrections for dust attenuation 
(Driver et al. 2008). The agreement is reasonable with the 
discs agreeing with the MGC data to within the error and 
the spheroid mass over-predicting the MGC value by a mod- 
est amount. It is worth noting from Fig. [5] that the stellar 
mass of spheroids is actually declining, with mass-loss hav- 
ing exceeding mass-gain for the past 9 billion years. For the 
discs, the two almost exactly balance, such that the over- 
all stellar mass density appears to asymptote to a constant 
value around the present day. 

Also shown in Fig. [S] as grey shaded data are the com- 
pendium of total stellar mass estimates given in Wilkins, 
Trentham & Hopkins (2008a). These data clearly fall signif- 
icantly below the black shaded line, highlighting a significant 
discrepancy between the total stellar mass inferred from the 
cosmic star-formation history and that derived from direct 
empirical constraints. This offset is well known and discussed 
in detailed in Wilkins, Trentham & Hopkins (2008a), here 
we make two additional comments: (1) the shape of the data 
and the black curve do broadly agree with a x 2 offset at al- 
most all ages, (2) the z = data from the MGC includes 
detailed dust corrections for both the optically thin and op- 
tically thick regions and is typically a factor of x2 higher 
than most local measurements. It is possible then that the 
values from the literature are missing mass embedded in op- 
tically thick regions. Perhaps a more likely explanation, also 
put forward by Wilkins et al. (2008) is that the IMF was 
simply lighter at earlier times. This would reconcile quite 
nicely as the low-z CSED is fairly impervious to the very 
low mass-end of the IMF. 



5.2 Dependency on assumed IMF 

We now briefly explore the impact of the adopted IMF. 
Fig. |4] shows the CSED predictions for the evolving metal- 
licity scenario using either a (top) Salpeter (1955), (centre) 
Kroupa et al. (1993), (bottom) or Kroupa (2001) IMF. Note 
that the Kroupa (2001) IMF is extremely close in form to 
the Chabrier (2003) IMF and therefore the Kroupa (2001) 
prediction can be taken as representative of both. 

Essentially all IMFs provide an equally good fit to the 
CSED except in the UV regime. At wavelengths longer than 
M-band (> 0.4/^m) the resulting shape is not particularly 
sensitive to the detailed shape of the IMF. This is mainly 
because at z = stars close to solar luminosity, where the 
IMF is least contentious, are dominating most of the CSED. 
Systems at the very high-mass end which formed at high 
redshift will no longer be contributing to the CSED whereas 
very low-mass stars are yet to dominate the near-IR flux. 
The CSED is therefore unable to constrain the IMF other 
than the normalisations required for these IMFs are gener- 
ally higher than for the models based on BG03. 



5.4 Discussion 

At this point we have a simple heuristic model which adopts 
two simple axioms motivated by the physically distinct ap- 
pearance of spheroids and discs in the nearby Universe. 
These axioms combined with the empirical compendium of 
AGN and cosmic star- formation activity/history are able to 
reproduce the z — CSEDs of spheroids and predict the 
mean mass and metallicity evolution of present day discs 
and spheroids. The model also provides a complete descrip- 
tion of the energy output from stars within those systems 
which will eventually make up the local spheroid population 
(projenitors) as a function of redshift, the metallicity build- 
up, and suggests key cross-over epoch at z ~ 1.6 between 
the hot and cold mode evolution. This later transition red- 
shift is consistent with the obvious change in morphologies 
seen in HST images at this redshift (e.g.. Driver et al., 1998, 
figure 3). 

However an obvious weakness is that the model provides 
no clear prediction of the morphological, size- and shape- 
evolution, merger rates, or the clustering of the galaxy pop- 
ulation. Furthermore the model does not stipulate the ac- 
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tual mechanism by which star-formation is occurring and for 
the hot mode could be some combination of monohthic col- 
lapse, major merging, and/or clump migration. Within the 
recent literature the exact status of the z ~ 2 population 
is also unclear. Deep IFU studies (e.g., Forster Schreiber et 
al., 2009,2011) find that the majority of star-formation at 
2^2 appears to be taking place in rotating clumpy disc 
structures with no obvious central bulge component. Simi- 
larly Chevance et al., (2012, see also Weinzirl et al., 2011) 
from a study of 31 high-z galaxies, find that the Sersic in- 
dices are significantly flatter than one would expect for low-z 
spheroids and obvious disc structures are present in many 
cases. From our Fig. 1 we can see that at 2; --^ 2 we are still 
within the epoch where spheroid formation should be domi- 
nating. However spheroid formation does not necessarily im- 
ply spheroid morphologies until after some unspecified time- 
lag in which the system settles. In fact violently star-forming 
systems will inevitably appear blue, asymmetrical, gas-rich, 
and dusty, i.e., quite disc-like in many aspects. Other stud- 
ies, e.g., van Dokkum (2008), find that 45% of massive 
galaxies at z ~ 2.3 do indeed have evolved stellar popula- 
tions, little or no ongoing star- formation, and compact early- 
type morphologies. Hence a picture of a spheroid population 
emerging from a highly turbulent progenitor phase around 
z ~ 2 appears to be qualitatively consistent with our model. 
Alternatively our model may need to be adjusted to allow 
for gas infall and disc-formation from the outset with some 
fraction of the disc- formed stars merging into bulges, i.e., a 
relaxation of the ma:ximal spheroid formation axiom. This 
would have the net effected of also increasing the cross-over 
redshift to > 1.6 

A further intriguing observations is that high-z 
spheroids are significantly more compact than nearby ellip- 
ticals by factors of x3 — 4 at fixed stellar mass (e.g., Daddi 
et al., 2005; van Dokkum et al., 2008 etc). Within our sce- 
nario this could be consistent with the high-z sample being 
"naked" -bulges yet to grow discs or yet to be "puffed-up" 
through successive minor merger interactions or adiabatic 
expansion. These two pathways, disc-growth verses "puff- 
ing", are likely to be strongly environmentally dependent 
with minor mergers more frequent in high density environ- 
ments, and gas infall more prevalent in low-density environ- 
ments. A particularly interesting comparison might there- 
fore be the mass-size relation of high-z spheroids to low-z 
bulges. 

With the caveat that the morphology and size evolution 
within our model is unspecified we nevertheless appear to 
have a prediction of the energy output of spheroid and discs 
projenitors over all epochs (Fig.|3]top panel), the mean gas- 
phase metallicity history for each population (Fig. [2]), and 
the build-up of stellar mass (Fig.[S}. Whether one can readily 
distinguish these populations observational however is an 
open question. 

Finally it is worthwhile reiterating that this model con- 
tains no tunable parameters nor any dependency on initial 
conditions beyond the underlying cosmology. The model is 
built entirely from empirical data and provides a fully con- 
sistent empirical scaffolding upon which more physically mo- 
tivated models can be built. Our conclusion is that the ini- 
tial axioms on which the model is based are viable and the 
star-formation histories defined are tenable. 

Further studies of the variation of the z = CSED and 




Time since Big Bang (Gyrs) 

Figure 6. As for Fig. [T] (lower) except with the spheroid star- 
formation history down-weighted by 25% by incorporating the 
CSED constraints from Fig. [3] (upper). 

its dependency on environment should enable an investiga- 
tion into dependencies on clustering, and to assess whether 
star-formation proceeds more rapidly or whether it is merely 
the relative mix of spheroid v disc formation which is chang- 
ing. Similarly, using observations at intermediate redshift, it 
should be possible to compare data from high-z studies to 
the predictions of our two-phase model. Both of these av- 
enues will be explored in future papers. 



6 CONCLUSIONS 

From two very simple axioms: (1) that AGN activity traces 
spheroid formation, and (2) that the CSFH is dominated by 
spheroid formation at high redshift, we are able to derive 
simple expressions for the cosmic star-formation histories of 
spheroids and discs. Following comparisons to the z = 
CSED for spheroids and discs we find a modest downward 
adjustment of 25% provides the optimal fit resulting in our 
final recommended star-formation histories of: 

PS = 50.77 X iQ-^hUJl^f-^eM-j^) (4) 
pn^il.m X lQ--^hU-^^r^°eM-T^^) (5) 

tGyrsAlO.7 IGyrsfto.7 

where ^ is the IMF multiplier as given in Table [1] Fig. [6] 
shows these final relations compared to the compendium of 
data provided by Hopkins & Beacom (2006) and despite 
the renomalisation of the spheroid star-formation history 
still provide a perfectly satisfactory description of the global 
CSFH. 

Adopting a Baldry & Glazebrook (2003) IMF and using 
these expressions to predict the z = CSED, we are able to 
provide a satisfactory explanation of the observed CSEDs of 
spheroids and discs from the FUV to the Jf-band. 
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The corollary of the simplicity of the two-phase model, 
however, is that it lacks any prediction of the clustering sig- 
nature, environmental dependencies, or the merger histories, 
although these can be built in at a later stage. Perhaps the 
key gain, in an era of hidden tunable parameters, is that with 
the adoption of a universal IMF and a stellar evolution code 
there are essentially no free parameters. Strictly speaking 
this is not precisely true as the detailed modelling of stellar 
evolution typically comes with options and there is arguably 
a choice of IMFs and also whether it is Universal or varies 
over cosmological time (see for example Wilkins et al. 2008a, 
2008b; Gunawardhana et al. 2011). On this last subject of 
the IMF it is worth reiterating that longward of 0.4/^m the 
z = CSED is not sensitive to the high-mass shape of the 
IMF (unless taken to the extreme). This is because at al- 
most all wavelengths, in the declining star-formation era to- 
day, the CSED is dominated by either the tip of the main 
sequence, which lies just below a solar mass and well above 
the mass range of contention (mid-optical to NIR), or the 
most recently formed stars (FUV to mid-optical). 

As a byproduct, the two-phase model also provides the 
CSED of spheroids and discs at every epoch in the Universe, 
along with the prediction of a clear-cut transition redshift 
at around z ~ 1.7 where galaxy evolution switches from 
evolution being dominated by major mergers to evolution 
being dominated by cold gas infall. Future work will include 
a broader wavelength baseline, bulge-disc decompositions, 
inclusion of the AGB energy output, and development of 
the model via comparisons to selected external data. 
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